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Abstract 

An aerodynamic sound generated by a flow inside a pipe is one of the noise problems. In addition to the aerodynamic sound 
sources, characteristics of the generated sound are sometimes affected by acoustic and vibration characteristics of the pipe. In 
this paper, we have performed experiments and simulations to clarify acoustic and flow induced sound characteristics of a 
T-shaped pipe with a square cross-section. The experiments and simulations were performed under several inflow velocity 
conditions. The results show that the characteristics of the flow induced sound in the pipe are strongly affected by the acoustic 
characteristics of the pipe interior sound field. The frequency characteristics of the sound sources generated by the flow depend 
on the inflow velocity conditions. However, that of the generated sound does not depend on the inflow velocity conditions, but 
on the acoustic frequency characteristics of the pipe. 
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Introduction 

An aerodynamic sound generated by a flow inside a pipe or a duct is one of the noise problems in the equipment 
such as air conditioner. An air conditioner consists of a duct, fan and openings. There are three types of flow in- 
duced sound sources in the air conditioner. One is the sound sources generated by the fan, so-called fan noise. The 
second is the sound sources generated by the flow around bending or expanding parts of the duct. The third is the 
sound sources generated by the flow near openings. However, it is still incompletely understood that which type 
of the sound sources is dominant. In order to reduce the noise, it is necessary to know the dominant sound sources. 

Studies on the aerodynamic sound generated inside the duct have been done experimentally. Watanabe et al. in- 
vestigated the acoustic power level generated by the air flow in duct systems [1]. Itamoto and Shiokawa studied 
flow noise and sound characteristics of glass fiver ducts [2], [3]. Their research results are also used as data collec- 
tion. Ishihara investigated acoustic and flow noise characteristics of straight and bending duct with some openings 
[4], [5], [6]. It is clarified that the noise reduction effect of the absorbing treatment could be neglected when noise 
sources were generated only by the flow. The noise reduction effects due to bending and use of inner guides or 
airfoils were clarified. In the case of pipe noise, Hambric et al. studied the flow-induced vibration and acoustic 
power spectra of a 90° piping elbow and proposed a procedure of the coupling of CFD and structural-acoustic 
models [7]. Zhang et al. studied the flow-induced vibration and noise of a 90° piping elbow and the effect of guide 
vanes [8]. They computed the hydrodynamic noise based on a hybrid LES/Lighthill's acoustic analogy [9] and clar- 
ified the effective position of guide vanes in reducing vibration and flow-induced noise in the 90° piping elbow 
with water. However, to the author's knowledge, there have not been many studies focused on effects of both 
acoustic and vibration characteristics of pipes on aerodynamic sounds. 

The purpose in this paper is to clarify the effect of the acoustic and vibration characteristics of the pipe on the aer- 
odynamic sound. We have performed both experiments and simulations of aerodynamic sounds generated by the 
flow in the T-shaped pipe with the square cross-section. To predict aerodynamic sounds computationally, the hy- 
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brid CFD/BEM method [10], [11] is applied. The mode shapes are calculated by means of FEM, and the acoustic 
characteristics are calculated by means of BEM including the structure-acoustic coupling effect. 


Numerical Procedure 


Transient CFD Simulation 


The transient flow fields in the T-shaped pipe with the square cross-section of a length D = 100mm on a side are 
simulated at Reynolds numbers. Re = 41075 and 78728. Inflow velocities are U = 6 and 11.5m/s. Mach numbers are 
M=0.0177 and M=0.0338, respectively. Figure 1 shows a computational domain and schematic diagram. For these 
simulations, the three-dimensional computational domain in the T-shaped pipe has been applied, as shown in Fig.l. 
Unsteady flow fields are calculated using the commercial CFD code ANSYS Fluent version 15.0 and its incom- 
pressible LES (Dynamic Smagorinsky model) calculation features. The origin of the Cartesian coordinate is placed 
at the midpoint of the pipe. The domain contains 544,000 hex cells and 573,221 nodes. The cell spacing adjacent to 
the wall is 0.8mm (0.008D). Steady velocities are imposed on the inflow boundary. Pressure boundary conditions 
are applied on the outflow boundaries. No-slip conditions are applied on the other walls. Steady state simulations 
were performed using Spalart-Allmaras (S-A) turbulence model and then used as initial conditions of transient LES 
simulations. The transient simulations were performed for 7,800 time steps with a time step size At = le-4s, which 
corresponds to a Courant-Fredrichs-Lewy (CFL) number of 2.1. 

Inflow 

_ 

X. 

Outflow 

FIG. 1 DOMAIN SHAPE AND SCHEMATIC DIAGRAM 



Lighthill Equation 


Lighthill equation [12], [13] in the frequency domain is derived from the equation of continuity and com- 
pressible Navier-Stokes equation and as follows: 


(V 2 -I- k 2 )p = - 


d 2 Ti m 

dxidx m 


(i) 


where p is the acoustic pressure, k is the wave number, c is the speed of sound, l and m indicate each direction 
in the Cartesian coordinates, v is the flow velocity. T im is the Lighthill stress tensor and as follows: 

Tim — P^l^m T (P ~ C p)^lm ~ T Im / (2) 


where p is the density and 1.225kg/m 3 , 8 \j is the Kronecker delta and T im is the viscous stress tensor. For a low 
Mach number and high Reynolds number flow regime, the second and third terms of equation (2) are negligible [8], 
[10], [14]. Therefore, the first term is used for the present work. 


Extraction of Acoustic Source 

To convert the acoustic source time histories into the frequency spectra, the discrete Fourier transform (DFT) has 
been applied. The acoustics sources are extracted from 512 steps (from t = 0.52s to 0.6222s) flow field data, the sam- 
pling period is 2e-4s. 


Acoustic Calculation 

The acoustic characteristics are solved using the commercial BEM code WAON. In this solver, the following simul- 
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taneous linear equation is solved: 


(£* + B + C)p — jcopAv + p d . 


(3) 


Here, p is the acoustic pressure vector, v is the particle velocity vector and the entries of the influence coefficient 
matrices are represented as follows: 

Etj=\s tj , (4) 

A tJ = f r . v Nj(r q )G(r t ,r q )dS q , (5) 

BiJ = l r (6) 

c tj = j;k A N j(rM r t’ r 'i) ds <i’ ( 7 ) 

where S t j is Kronecker delta, and T is the total boundary. T v is a part of T where the surface is assumed to be 
vibrating. T A is also a part of T where the surface absorbs acoustic wave. r t is a position vector of the i -th node, 
r q is a position vector of the source point q, and Nj is the interpolation function of the ;-th node. 

With the number of nodes N, the component p of the vector p is as follows: 

p(r q ) = 11^=1 Nj(r q )pj. ( 8 ) 


The component p d of the vector p d is the direct pressure contribution from the acoustic source, which is evaluated 
by the following equation: 


Pd(r P ) 


1 9 2 roo ■T; m (r s , & ))e^l r P-^l ^ 

4ndxidx m J-°° |r p -r s | 1 


(9) 


where d 2 / dx x dx m is the directional derivative and V is the volume of the flow field (in this case, the region filled 
by CFD cell). r p is a position vector of the monitor point p, and r s is a position vector of the source point s. 


There are 2,780 boundary elements. The acoustic sources are extracted from CFD results, whose numbers are 
equivalent to the number of grids of the CFD model. Figure 2 shows the boundary elements. To concern acoustic 
waves moving from the inflow boundary to the outside, an impedance boundary condition is imposed at the in- 
flow boundary. The value of the impedance is pc. In the structure-acoustic coupled model, all surfaces other than 
inflow boundary are boundaries where the structure-acoustic coupling effect is considered. In the uncoupled mod- 
el, they are assumed to be rigid. 


Absorption boundary 
(Inflow boundary) 


* 

L. 

FIG. 2 BOUNDARY ELEMENT MESH 


Structural-Acoustic Coupling 

In this paper, the structural-acoustic coupled model is solved. In order to solve structural-acoustic coupled 
model, equations corresponding to the acoustic field and structure must be solved simultaneously. 

Displacement of the structure is described using equation of motion as follows: 

(—oo 2 M + jcoC + K)u = f s , (10) 

where M is the Mass matrix, C is the Damping matrix, K is the Stiffness matrix and f s is the excitation force 
vector. 
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With modal coordinates. Equation (10) is expressed as follows: 

( — 60 2 + 2jho)O) 0 + 0)1) q = (p T f s , (11) 

where h is the modal damping ratio, O) 0 is the matrix in which eigenvalues are arranged diagonally, <p T the ei- 
genvector transposition, and q the modal participation factor. 

In order to consider the acoustic contribution to the structural model, an additional force term from the acoustic 
field is introduced as follows: 


(- 0 ) 2 + 2jho)O) 0 + <a%)q = <p T (j s - ps ), (12) 

where s is the area affected by the acoustic pressure. 

On the other hand, to consider the structural contribution to the acoustic model, an additional velocity term from 
the structure is introduced to Equation (3): 

(£* + B + C)p = jcopAv + jcopA v s + p d , (13) 

here v s the vibration velocity vector of the structural model and. 


A'ij = fr 3 N j( r M r i- r ‘i)dS q , (14) 

where r 3 is the boundary on which structural-acoustic coupling effect is considered. T 3 corresponds to all of the 
surfaces other than the inflow boundary in Fig.2. 


The term v s in Equation (13) is modified into an expression in which modal coordinates are used. Combining 
Equation (12) and (13), we obtain the following coupled equation of acoustic and structural models: 


E + B + C o) 2 pA (p rP| _ (jcopAv + p d j 

-(p T s -a) 2 + 2jho)O) 0 + o>q _W l <p T f s ) 


(15) 


The modal information of the pipe has been extracted using ANSYS Mechanical version 15.0. Figure 3 shows 
boundary conditions in the modal analysis. 



FIG. 3 BOUNDARY CONDITIONS IN THE MODAL ANALYSIS 

Young's Modulus, Poisson's ratio and density of the T-rectangular pipe are respectively, 3.14GPa, 0.35 and 
1190kg/m 3 . The damping ratio is 0.05 (5%). The thickness of the pipe is 5mm. As the result of the modal 
analysis, typical mode shapes are obtained, as shown in Fig. 4. 



10th mode at /= 190.31Hz 11th mode at /= 190.35Hz 
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21st mode at /= 353.76Hz 22nd mode at/= 353.76Hz 

FIG. 4 TYPICAL SHAPE MODES OF MODAL ANALYSIS 





Experimental Setup 

Experimental setup is shown in Fig. 5. With a microphone that is placed at 100mm from the centre of the outflow 
boundary, sound pressure signals are measured and converted to frequency domain with FFT Analyser. The fre- 
quency range is from 20Hz to 1000Hz, the sampling frequency is 10000Hz, the number of averages is 1000, and the 
frequency resolution is 10Hz. Inflow velocities are about 6, 8, 10, 11.5m/s. Background sound pressure level is be- 
low 30dB. 



100mm 

FIG. 5 EXPERIMENTAL SETUP 


Results and Discussion 

Transient CFD Results 

Instantaneous snapshots of the vorticity field at Z = 0 plane are shown in Fig. 6. Vortices generated near corners of 
the pipe are shed from the corner and travels downstream. Shed vortices merge and interact with the pipe wall. 
Motions of these vortices play an important role of the sound generation. Figure 7 shows instantaneous snapshots 
of the Lighthill stress tensor pv l v m at Z = 0 plane. Significant acoustic sources exist near corners and the wall of the 

pipe in Fig. 7, as in the vorticity field. Figure 8 shows the magnitude of the volumetric acoustic sources (Lighthill 
stress tensor) extracted from the CFD results. As shown in Fig. 8, magnitudes of the acoustic sources are relatively 
symmetry with respect to the YZ plane. The frequency spectra of the volumetric Lighthill stress tensor pv l v l 
at (ID, 0, 0) are shown in Fig. 9. Figure 9 (a) shows the spectra of normalized frequency (Strouhal number, 
St = fD/U), and Figure 9 (b) shows the spectra of non-normalized frequency. Some outstanding peaks are sur- 
rounded by blue and red circles, as shown in Fig.9 (a) and (b). The Strouhal numbers of the peak surrounded by 
the blue and red circles are St=1.78 and 3.74, respectively in Fig.9 (a). The Strouhal number, St=1.78 corresponds 
with 205.1 [Hz] at U=11.5m/s and 107[Hz] at U=6m/s in Fig.9 (b), respectively. In the same way, the Strouhal num- 
ber, St=3.74 corresponds with 429.9[Hz] at U=11.5m/s and 224.3[Hz] at U=6m/s, respectively in Fig.9 (b). Therefore, 
it is clear that the peak frequencies of the acoustic sources depend on the inflow velocity. 
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Time = 0.56560 






FIG. 6 INSTANTANEOUS SNAPSHOTS OF VORTICITY FIELD AT Z=0 PLANE (LEFT: 11.5 m/s, RIGHT: 6 m/s) 



0 0006+000 


D 



FIG. 7 INSTANTANEOUS SNAPSHOTS OF LIGHTHILL STRESS TENSOR pv t v m AT Z=0 PLANE (LEFT: 11.5 m/s, RIGHT: 6 m/s) 
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FIG. 8 MAGNITUDE OF VOLUMETRIC ACOUSTIC SOURCES pV t V m AT Z=0 PLANE (LEFT: 11.5 m/s, RIGHT: 6 m/s) 
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(a) NORMALIZED FREQUENCY (STROUHAL NUMBER) 



(b) NON-NORMALIZED FREQUENCY 

FIG. 9 SPECTRUM OF LIGHTHILL STRESS TENSOR pv^ PER UNIT VOLUME 


Far-Field Acoustics 

To investigate the effect of the structural-acoustic coupling on the far-field sound pressure, we have calculated 
sound pressure fields in both cases of the coupled and uncoupled model using the acoustic sources extracted from 
the CFD results. Figure 10 shows the comparison of sound pressure levels in the cases of the coupled and uncou- 
pled model. The position of the monitor point is shown in Fig.5. Peak levels of the sound pressure in the case of the 
coupled model are smaller than those in the case of the uncoupled model. However, there is no significant differ- 
ence between the coupled and uncoupled model in terms of the frequency characteristics of the far-field sound 
pressure. It indicated that there is relatively small effect of the structure vibration on the far-field sound pressure in 
the present work. 



FIG. 10 COMPARISON OF SOUND PRESSURE LEVELS IN THE CASES OF THE COUPLED AND UNCOUPLED AT U = 11.5m/s. 
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Figures 11 and 12 show the frequency spectra of SPL at the position. Both the simulated and measured data are 
shown in Fig. 11, and only the measured data are shown in Fig. 12. Comparing the simulated and measured spec- 
tra, the simulated spectra have some intense peaks. It is likely that the difference between the measured and simu- 
lated is attributed to the acoustic damping. Therefore, we have calculated the structural - acoustic coupled re- 
sponses with the weak sound absorption on all boundaries except the inlet. The absorption coefficient is 0.05 in the 
present work. Figure 13 shows the frequency spectra of SPL with the weak sound absorption. The intense peaks are 
suppressed. There is a possibility that an acoustic damping in addition to the structural damping acts in the ex- 
periment. In regard to the acoustic damping, it is necessary to consider for our future work. 

As shown in Fig.9, the frequency characteristics of the acoustic sources depend on the inflow velocity. However, 
both the measured and simulated frequency characteristics of SPL do not depend on the inflow velocity, and peak 
frequencies of SPL are almost the same between 11.5 m/s and 6 m/s cases, as shown in Figs. 11, 12 and 13. Figure 14 
shows contours of the sound pressure field. Depending on the inflow velocity, magnitudes of the sound pressure 
are different between the two cases. However, directivities of the sound pressure do not depend on the inflow ve- 
locity, and they are almost the same in the two cases. Table 1 shows the overall sound pressure level (OASPL) 
(50-1000Hz). The OASPLs of the measured and simulated data show a relatively good agreement. 





FIG. 13 SPECTRUM OF SPL NEAR THE OUTFLOW BOUNDARY WITH WEAK ACOUSTIC ABSORPTION 
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FIG. 14 SOUND PRESSURE FIELD AT Z=0 PLANE (LEFT: U = 11.5m/s, RIGHT: U=6m/s) 


TABLE 1 OVERALL SOUND PRESSURE LEVEL (50HZ-1000HZ) 


Inflow Velocity 

Measured 

Simulated 

Simulated (Weak Absorption) 

6m/s 

90.1dB 

91.8dB 

89.22dB 

8m/s 

95.3dB 

— 

— 

lOm/s 

100.2dB 

— 

— 

11.5m/s 

104.3dB 

109.8dB 

104.96dB 


Acoustic Characteristic of the Pipe 

As shown in Fig. 12, peak frequencies of SPL do not depend on the inflow velocity. Therefore, it can be presumed 
that the acoustic modes of the pipe interior sound field strongly affect the frequency characteristics of the flow in- 
duced sound in the pipe. To clarify the acoustic characteristic of the pipe, the acoustic frequency responses have 
been calculated using the monopole point sources (without the flow) shown in Fig.15. To examine the position de- 
pendence of the sound source, two types of analysis models are used, as shown in Fig.15. In the case of Type 1, the 
point source is located at the origin ("Symmetric source" position case). In the case of Type 2, the point source is 
located at (4D, 0, 0) ("Asymmetric source" position case). The magnitude of the point sources is IPa in all frequen- 
cies. Figure 16 shows acoustic frequency charecteritics of the pipe (without the flow). The measured data (with the 
flow) is also plotted in Fig.16. Some outstanding peaks are surrounded by red and blue circles in Fig.16. Some 
acoustic modes of the pipe are also shown in Fig. 17. Figure 16 and Figure 17 show that in the case of Type 2, more 
acoustic modes are excited than in the case of Type 1. Comparing the cases of Type 1 and Type 2, in the case of 
Type 1, the simulated data is smooth and similar to the measured data (with the flow). However, in the case of 
Type 2, some acoustic modes (peaks surrounded by the red circle) appear and their modes also appear in the 
measured data. The results show that the acoustic frequency characteristics of the pipe depend on the position of 
the acoustic source, and are similar to the frequency characteristics of the flow induced sound in the pipe in terms 
of the peak frequency. Also, the results explain that the frequency characteristics of the flow induced sound in the 
pipe are affected by the acoustic frequency charecteristics of the pipe. 
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(a) Typel ("Symmetric source" case) 
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(b) Type2 ("Asymmetric source" case) 


FIG. 15 ACOUSTIC SOURCES FOR THE FREQUENCY RESPONSE ANALYSIS 



FIG. 16 ACOUSTIC FREQUENCY CHARACTERISTIC OF THE PIPE (WITHOUT THE FLOW) 
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FIG. 17 ACOUSTIC MODE OF THE PIPE (LEFT: TYPE1 (SYMMETRIC SOURCE), RIGHT: TYPE2 (ASYMMETRIC SOURCE)) 


Conclusions 

In this paper, we have performed the experiments and simulations on the acoustic, vibration and aerodynamic 
sound characteristics of the T-shaped rectangular cross-sectional pipe. The experiments and simulations were per- 
formed under several inflow velocity conditions. The following concluding remarks are obtained in this paper. 

(1) To investigate the effect of the structural-acoustic coupling on the far-field sound pressure, the struc- 
ture-acoustic coupled and uncoupled model are calculated using the acoustic sources extracted from the 
CFD results. The effect of the structure vibration on the far-field sound pressure is relatively small in the 
present work. 
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(2) Magnitudes of the generated aerodynamic sound pressure vary depending on the inflow velocity. However, 
the directivities and frequency characteristics of the sound pressure do not depend on the inflow velocity, 
and they are almost the same for all inflow velocities. 

(3) The acoustic frequency responses of the pipe interior sound field were calculated using two types of the 
monopole point sources (without the flow). The acoustic frequency characteristics of the pipe depend on the 
position of the acoustic source. In terms of the peak frequency, they are similar to the frequency characteris- 
tics of the flow induced sound in the pipe. The results explain that the frequency characteristics of the aer- 
odynamic sound in the pipe are affected by the acoustic frequency charecteristic of the pipe. 

The further study of the case that the thickness of the pipe is reduced is underway to clarify the effect of the struc- 
tural-vibration on the far-field sound. 
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